Shear- Transformation-Zone Theory of Glassy Diffusion, Stretched Exponentials, and 

the Stokes-Einstein Relation 
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The success of the shear-transformation-zone (STZ) theory in accounting for broadly peaked, 
frequency-dependent, glassy viscoelastic response functions is based on the theory's first-principles 
prediction of a wide range of internal STZ transition rates. Here, I propose that the STZ's are 
the dynamic heterogeneities frequently invoked to explain Stokes-Einstein violations and stretched- 
exponential relaxation in glass-forming materials. I find that, to be consistent with observations of 
Fickian diffusion near T g , an STZ-based diffusion theory must include cascades of correlated events, 
but that the temperature dependence of the Stokes-Einstein ratio is determined by an STZ-induced 
enhancement of the viscosity. Stretched-exponential relaxation of density fluctuations emerges from 
the same distribution of STZ transition rates that predicts the viscoelastic behavior. 



Among the deepest challenges in glass physics is un- 
derstanding the temperature dependence of the Stokes- 
Einstein ratio in the neighborhood of the glass temper- 
ature T 9 .[ll-Q An apparently related phenomenon is the 
stretched-exponential relaxation of density fluctuations 
and other correlations. Both phenomena have been cited 
as evidence for dynamic heterogeneities in glassy systems, 
which supposedly provide rapid diffusion paths relevant 
to the first case, and a variety of localized, relaxation en- 
vironments in the second. [3] Here, I propose that these 
heterogeneities are the same as the shear-transformation 
zones (STZ's) that Bouchbinder and I Q invoked to 
understand the broad range of time scales observed in 
frequency-dependent viscoelastic response functions. @ 

As described in earlier publications STZ's are 

localized, shearable, thermally activated, structural fluc- 
tuations that appear and disappear on time scales of the 
order of the a relaxation time T a . In the presence of 
an applied shear stress, they are the places where driven 
shear rearrangements occur. They must also be associ- 
ated with the self diffusion of a tagged molecule. At most 
times, such a molecule remains in its cage, surrounded 
by its neighbors. It participates in a rearrangement in- 
volving these neighbors only when a thermally activated 
structural fluctuation - a shear transformation - occurs 
at its location. When that happens, the tagged molecule 
jumps some distance, perhaps an intermolecular spacing 
a, in some direction, while one or more of its neighbors 
jumps in the other direction. A sequence of such jumps 
becomes a diffusive trajectory. 

For present purposes, consider only systems in thermo- 
dynamic equilibrium near their glass temperatures, and 
their linear responses to small deviatoric shear stresses s. 
The basic physics of the STZ model in this limit is ex- 
pressed by the master equation for the number densities 
n±{v) of STZ's oriented parallel or antiparallel to s: 



h±(v) = R(±s)n^(v) — R(=fs)n±(v) 
+ p(0)[(l/2)e- B *' e -n±(v)]. 



Here, the R(±s) are the rates of internal STZ transitions, 
9 = ks T is the temperature in energy units, and ez is the 
STZ formation energy. p{9) is the super- Arrhenius rate 
factor that rapidly decreases toward zero as 9 decreases 
through the glass temperature 9 g . Times are measured 
in microscopic units, perhaps molecular vibration peri- 
ods. The symbol v denotes the internal transition rate, 
v = 2 R(0); and writing Eq.([T]) with this explicit v de- 
pendence implies that a range of different values of v will 
be relevant to the phenomena of interest. Denote the 
probability distribution for the v's by p{v). 

To start, compute the steady-state plastic strain rate 
-y pl by setting h±{v) = in Eq.(Q}. The result can be 
written in the form 



where 
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There are several aspects of these equations that need 
emphasis. First, they assume that eo, the volume of the 
plastic core of a shear transformation in units of a 3 , is 
independent of v. Second, Eq.(|3]) includes a definition 
of r Q that is not tied directly to the temperature depen- 
dence of the viscosity r\\ that is, the 9 dependence of rj 
will not be simply proportional to r a {9). Third, T(s) 
is the stress-induced bias between forward and backward 
transitions. The fact that this quantity is equal to vq s/6 
(where vo ~ a 3 ) is the classic Stokes-Einstein formula. 
In [10(, this result was shown more generally to follow 
from the second law of thermodynamics, with 9 equal to 
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the thermodynamically defined effective temperature \- 
A simple application of the preceding analysis to sheared 
systems, where the transverse diffusion constant is a 2 /r Q , 
confirms that x is the same as the effective temperature 
determined by a fluctuation-dissipation relation, at least 
within the STZ theory. (See and references cited 
there.) Throughout the following, thermodynamic equi- 
librium implies that x — 

In [5(, Bouchbinder and I wrote the ^-distribution in 
the form 



p(v) 



A 



v {{yjv*)^ + {u*/v)^} 



(5) 



where A is a normalization constant. The large- f behav- 
corresponds to a Boltzmann distribution 



lor, p • 



over barrier heights, roughly consistent with that found 
experimentally by Argon and Kuo.[l2j We used £1 ~ 1, 
and C ~ 0.4 in accord with viscoelastic data for various 
materials near their glass temperatures. The cutoff rate 
v* was determined by noting that the internal transition 
rate v cannot be slower than the rate of spontaneous 
structural rearrangements. Therefore, v* ~ t^ 1 , which, 
for systems near their glass temperatures, is a factor of 
about 10 3 smaller than p. These first-principles estimates 
of p{y) and v* were accurately confirmed by the experi- 
mental data. 

Equation ([3]) implies that there are two different lim- 
iting populations of STZ's. There are "fast" STZ's, with 
v p, for which the time scale t{v) « r a . These STZ's 
may make multiple, back and forth transitions between 
their two states during their lifetimes. There is also a 
substantial population of "slow" STZ's with v <C p and 
t{v) rs (p/v) r a 3> r a , which make at most one transition 
before disappearing. Both populations are important for 
determining the inverse viscosity: 
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pi 



dvp(v) 



(6) 



Turn now to the statistical problem of self diffusion. 
Here, instead of averaging over probabilities of single 
events as above, we must consider sequences of mul- 
tiple hopping events. The statistical weight assigned 
to any such sequence depends on how far the molecule 
moves and how long it takes to get there. Technically, 
the molecule executes a "continuous-time random walk" 
(CTRW). 0, H| Note that the model of STZ-enabled 
diffusion steps contains no information about dynamics 
on sub-molecular space or time scales. In effect, it is 
coarse grained on the scale of the molecular spacing, and 
it "knows" nothing about fast vibrational motions other 
than that they carry the thermal noise that activates 
barrier-crossing transitions. Therefore, this model can- 
not describe j3 relaxation or any details of how a molecule 
moves within the cage formed by its neighbors, or of how 
it escapes from that cage. 



The basic ingredient of the CTRW analysis is the one- 
step probability that, after waiting a time t, an STZ 
fluctuates into existence at the position of the tagged 
molecule, and the molecule jumps a distance z. (There 
is no loss of generality in projecting the diffusion mecha- 
nism onto one dimension.) For an STZ with an internal 
rate v, this normalized one-step probability distribution 
is 



1 



t{v) 



-t/T(» 



(7) 



where f(z, v) is the probability of a jump of length z. As 
indicated, f(z, v) may depend on v. The crucial assump- 
tion here is that the time t(v) appearing in this formula 
is the same as the quantity defined in Eq. ([3]) . In both 
cases, t{v)~ x is the average rate at which STZ's appear, 
make one net internal transformation (perhaps after mul- 
tiple back-and-forth transformations), and then disap- 
pear. This elementary fluctuation mechanism should be 
common to both plastic deformation and self diffusion. 

The probability of any random walk consisting of a 
sequence of such steps is a product of space-time convo- 
lutions of factors ip(t, z, v). If the STZ events are uncor- 
rected - a major assumption - then each factor can be 
averaged independently over v. The walk ends after a 
final time interval t within which the molecule does not 
make a further jump; thus the complete path probability 
contains one factor 



r(y) 



df e -*'M»0 = e -*AM ) (g) 



again included in the convolution integrals and averaged 
over v. 

These space and time convolutions are converted into 
products by computing Fourier and Laplace transforms, 
in terms of which the multi-step walk probabilities can 
be summed to all orders. The self-intermediate scattering 
function F(k,t), i.e. the fc'th Fourier component of the 
diffusion profile, is given by an inverse Laplace transform: 



F(k,t) 
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dw e wt/Ta K(w) 
2ni 1 - j(k,w) ' 



where 



dv m . , A(,) 



S > 0, (9) 
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o w + X(v) t(v) v + p 

is the Laplace transform of the averaged <f>(b, u), and 



J(k, w) 



dv p(y 



w + \{y) 



(11) 



is the analogous transformation of ij)(t, z, v). 

A crucial ingredient in these formulas is the jump- 
length distribution f(k,v). If, according to the discus- 
sion following Eq.Q, each STZ event moves the tagged 
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FIG. 1: Stokes-Einstein ratio Cse for ( = 0.4, 0.5, 0.6, 0.7, 
and 0.8 from bottom to top. 



molecule a mean-square distance a 2 , then it would be 
natural to choose a Gaussian distribution f(k,u) — 
exp(— k 2 a 2 /2) independent of v. With this choice, how- 
ever, the slow STZ's play the role of deep traps [l4|, ef- 
fectively shutting down long-time diffusion. To see this, 
compute the mean-square displacement z 2 (t) 7 which, for 
a z^-independent f(k), becomes 



d 2 F(k,t) 
dk 2 



dw e w */ r ° 



2ni 



w 2 K(w) 



(12) 



If p{v) cuts off sharply at small v : then K(w) is an- 
alytic at w = 0, and the long-time behavior is ob- 
tained by integrating around the pole there. The re- 
sult is that z 2 (t) ps V t, with V Q = {a 2 /r a ) A'(0) _1 . 
Therefore, ordinary diffusion is suppressed by a factor 
KiO)- 1 ~ v*jp < 1. For Ci < 1 in Eq.@, K{0) diverges, 
and T>q vanishes. This result is manifestly inconsistent 
with the experimental data reported in [l|-|3j , where Fick- 
ian diffusion was observed near the glass temperature. 

Given the p(v) determined by viscoelasticity @, this 
discrepancy between diffusion theory and experiment 
falsifies the postulate of uncorrelated short jumps. It 
seems, instead, to indicate that the cascades of correlated 
STZ "flips" seen in low-temperature numerical simula- 
tions (see [HI and earlier papers cited there) are rele- 
vant to laboratory systems at temperatures near T g . A 
systematic generalization of the CTRW analysis to in- 
clude strongly correlated sequences of jumps is beyond 
the scope of this investigation; but it is plausible that 
such sequences might appear approximately in the form 
of anomalously long jumps initiated by individual slow 
events. Unlike the fast STZ's, the slow ones are stable 
only in near-equilibrium glass formers at low tempera- 
tures, and they only occasionally undergo shear transfor- 
mations during their lifetimes. Their small values of v im- 
ply that their internal energy barriers are high; therefore, 
when they do make transitions, they release substantial 
amounts of energy that could trigger subsequent events. 



Suppose, as a first guess, that a slow event sends a tagged 
molecule on a random walk that lasts for a time t(i>), 
during which the jump rate is t" 1 . The mean-square dis- 
placement during this walk would be a 2 t(v) /t q , which 
is equal to a 2 for fast STZ's, but becomes substantially 
larger for slow ones. 

The corresponding jump-length distribution is 



f(k,u) 



-fc"V/2A(i/) _ & -k' 



r{u)/2; 



(13) 



This choice of f(k,v) solves the problem of overly slow 
diffusion. When it is used in Eq.©, the quantity K(w) 
in the denominator of Eq. (fT2| cancels out, and z 2 (t) = 
a 2 t/T a at all times. This means that diffusion in this 
approximation is normal in the sense that T> = a 2 jr a . 

With this expression for V, and Eq.©, the Stokes- 
Einstein ratio becomes 



a T>r] 



dvp(y) 



C S b{9), (14) 



where ao = eovo/a 2 . Since r a has cancelled out of 
Eq. (ll4[) . and p(y) is a function of v/v* , the only temper- 
ature dependence of Cse(9) occurs via the ratio v* / p ~ 
exp(— ez/0). In [B[, we estimated that exp(— ez/0 g ) ~ 
10~ 3 . That estimate has been used to plot the graphs 
of Cse as functions of T/T g shown in Fig. [T] The bal- 
ance between the contributions of slow and fast STZ's in 
Eq. (ri4"|) is determined largely by the exponent C in Eq.©, 
thus Cs e is shown in the figure for five different values of 
£. Note that this theory does not (yet) have enough mi- 
croscopic physical content to make a connection between 
low and high-temperature behaviors, or even to make 
material-specific predictions of the shape of p{v) or the 
length scale ao- We know that, at higher temperatures or 
in noisier driven situations, p(v) becomes sharply peaked 
at much higher values of v, so that Cse — 1- Thus, 
the values of CsEifig) in Fig. [Dimply Stokes-Einstein 
violations of roughly the magnitude seen experimentally, 
arising from the contributions of the slow STZ's to the 
value of the viscosity according to Eq.©. 

Finally, consider the self-intermediate scattering func- 
tion F(k,t) given by Eas.piH]). To evaluate F(k,t) 
numerically, close the contour of integration in Eq.© 
around the branch cut on the negative real w axis. 
Some results are shown in FigJ5] The parameters cor- 
respond roughly to a metallic glass near T g : p/v* = 
10 3 , C = 0.4, and Ci — 1- The wavenumbers are 
ka = 2.0, 1.0, 0.5, and 0.1, from left to right. All four 
curves exhibit stretched-exponential relaxation of the 
form exp[-c(t/r Q ) b ], with b ^ 0.43, 0.45, 0.58, and 0.78, 
in the same order of decreasing k. 

Several features of these results are notable. First, 
there is a Fickian limit at small k, where b approaches 
unity, and F(k,t) ~ exp (— T> k 2 1/2) as shown by the 
dashed line in Fig J5] for ka = 0.01. Second, there is a 
non-Fickian limit at large k, where the curves begin to lie 
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FIG. 2: Computed self-intermediate scattering functions for 
ka = 2.0, 1.0, 0.5, and 0.1, from left to right. The correspond- 
ing stretched-exponential indices are 6 = 0.43, 0.45, 0.58, and 
0.78. The dashed curve is the Fickian limit for ka = 0.01. 
Only a relaxation appears here. Space and time scales rele- 
vant to /3 relaxation have been averaged out. 



on top of each other, and b approaches 0.4. This limiting 
behavior occurs because any choice of the factor /(fc, v) 
vanishes at large fc, leaving the initial, fc-independent 
function K(w) in the numerator of Eq.® as the only 
contribution to the relaxation function. Note that 



K(t) = 



dw 
2ttz 



? wt / T «K(w) 



dv p(v) 



-t/rO) 



(15) 

has the form of a local relaxation function averaged over 
a distribution of environments with different relaxation 
times r{y). The transition between large-fc and small-fc 
behavior is temperature dependent via the quantity X(w) 
in Eq.©, which produces a 9 dependence here in much 
the way it did for the viscosity in Eq.©. At fixed fc, 
according to these equations, the system changes from 
stretched-exponential to Fickian as the temperature in- 
creases. It is not a coincidence that, in the large-fc limit 
given in Eq. (fl~5"j) and at low temperatures, b = £; but 
b — ( is not even an exact result for all choices of £ and 
£i, nor is it exact in the limit of long times. 

The analysis presented here calls for a unified investiga- 
tion of Stokes-Einstein violations, stretched-exponential 
relaxation, and frequency dependent viscoelastic re- 
sponse functions, all in comparable glass-forming liquids 
near their glass temperatures. There are major uncer- 
tainties, the most serious being the assumed /(fc, v) in 
Eq. (rr3l) . Without some such assumption, the existence 
of the slow STZ's deduced from the viscoelastic experi- 
ments would be strongly inconsistent with the observa- 
tions of Fickian diffusion near T g . [3-[3| However, it is not 
clear that Eq. (fT5|) is a realistic approximation, or whether 
a more accurate description of correlated diffusion steps 



might modify the theoretical Stokes-Einstein ratio. 

More generally, checking the consistency of the various 
formulas shown here would test the theory as a whole 
and its interpretation of dynamic heterogeneities. For 
example, it ought to be possible to measure viscoelastic 
responses (or internal-friction functions) using the same 
materials, under the same conditions, that are used for 
measuring F(k,t) in the large-fc limit shown in Eg. (1151) . 
Then, it should be possible either to back out consistent 
approximations for p(v), or else to falsify the predicted 
consistency. Similarly, it would be useful to make the 
latter kinds of measurements for the same systems in 
which Fickian diffusion is observed, and thus to test the 
approximation for f(k, v) in Eq.flTS]). Detailed studies of 
the k and 9 dependences of F(k,t), again on the same 
systems for which other measurements are being made, 
might further resolve the outstanding issues. 
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